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Abstract 

We define phylogenetic projective toric model of a trivalent graph as a generalization 
of a binary symmetric model of a trivalent phylogenetic tree. Generators of the pro- 
jective coordinate ring of the models of graphs with one cycle are explicitly described. 
The phylogenetic models of graphs with the same topological invariants are deforma- 
tion equivalent and share the same Hilbert function. We also provide an algorithm to 
compute the Hilbert function. 
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Introduction and background 



The inspiration for this work are toric varieties arising in computational biology, or more precisely in phylo- 
genetic algebraic geometry. The references to the subject include |Fel04] . |PS05] and |SS03) . 

Markov models on phylogenetic trees are statistical models describing evolution. They are usually defined 
as a subset of the probability simplex, parametrized by a subset of matrices depending on the model. Among 
them there are group-based models on phylogenetic trees. These are special, as their projective versions, that 
is the Zariski closure of the parametrization in the complex projective space, are projective toric varieties. 

Wc are interested in the simplest group-based models — binary symmetric models, also called the Jukes- 
Cantor models, on trivalent trees. The object of study is the generalization of those models to trivalent 
graphs. 

The article is organized in the following way: in Section 11.11 we give motivation to the subject and 
point out references where the main object of our study appears. Then we give a brief introduction to our 
main tools: we recall geometric invariant theory in Section 11.31 next, in Section 11.41 we set the notation for 
projective toric varieties and in Section 11.51 provide the description of GIT quotient of a projective toric 
variety by subtorus of its big torus. 

In Section [2. II we state combinatorial relations between the topological invariants of a trivalent graph and 
we prove that graphs with the same invariants are mutation-equivalent. In Section [2] we define the model of 
a trivalent graph as a GIT quotient of a product of P'^'s indexed by inner vertices of the graph. This implies 
that the model is toric and comes with an embedding into a weighted projective space. Our first result is 
Theorem 12.331 that lists the set of minimal generators of the projective coordinate ring of the model, when 
the underlying graph has the first Betti number at most one. Section |3] contains our second result — models 
of mutation-equivalent graphs are deformation equivalent. 

Theorem 13.51 Geometric models of connected trivalent graphs with n leaves and the first Betti number 
g are deformation equivalent in the projective toric variety Pg,n; which is a quotient of P " by a 

g-dimensional torus. 

In Section |4] we prove that the Hilbert functions of mutation-equivalent models are equal (Theorem l4.5|) 
and finally we compute these Hilbert functions explicitly. 

1.1 Motivation — Markov models on phylogenetic trees 

A phylogenetic tree is an acyclic connected graph with additional data attached to its edges and vertices. 
At a vertex v there is a finite ordered set Ay called an alphabet. At an edge with ends v and w there is a 
doubly-stochastic matrix (all rows and columns sums are 1) with the {i^jY^ entry indicating the probability 
of the i*'' letter of Ay being changed to the letter A^]- To construct a Markov model on a phylogenetic 
tree we first need to indicate a set of observable vertices, for example the leaves of the tree. Then the model 
is the subvariety of the probability simplex, parametrized by a subset of matrices that we only allow, given 
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by probabilities of observing letters at the observable vertices. We consider symmetric models, which means 
we allow symmetric matrices. Typically the observable vertices are the leaves of the tree. 

Apart form this real variety, one can consider its complex algebraic relaxation. That is, the parameters 
are allowed to vary in a complex projective space and we take the Zariski closure of the image. Then the 
model becomes a complex projective variety and can be studied by means of algebraic geometry. Binary 
symmetric models have additional structure — they are equipped with an action of a torus of dimension 
equal to the dimension of the model and thus they are projective toric varieties. This is an especially nice 
class of varieties, which have a combinatorial description by lattice polytopes. The geometry of the simplest 
group-based models — binary symmetric models with the restriction that the underlying tree is trivalent was 
the object of study of |BW07| . In that paper we described the corresponding lattice polytope and interpreted 
the models as a certain quotient of a product of three-dimensional projective spaces. 

1.2 Toric algebras of our graph models in the literature 

We generalize the quotient description of the tree models introduced in |BW07j and we again have a toric 
projective model, which this time is embedded in a weighted projective space. Such an embedding is always 
given by a graded lattice cone. We denote the cone for a graph Q by t{Q). 

The way we associate a lattice cone to a trivalent graph appears also in the work of Manon 'Man09! . He 
constructs a sheaf of algebras over the moduli stack Adg^n of genus g curves with n marked points and our 
semigroup algebras C[r(fJ)] are obtained by some initial term deformations from algebras above the most 
special points of A^g,™ in Manon's construction. 

Another place where our cones t{Q) appear is the Jeffrey and Weitsmann's |JW92) study of flat SU{2)- 
connections on a genus g Riemann surface. In their context the trivalent graph Q describes the geometry 
of the compact surface of genus g and thus has no leaves. A subset of Z-labellings of the graph, which 
are exactly points of our cone t(Q), are in 1-1 correspondence with the number of Bohr-Sommerfeld fibers 
which is the central object of study in |JW92| . By the Verlinde formula, the number of those fibers equals 
the dimension of holomorphic sections of powers of a natural line bundle on the moduli space of flat SU{2) 
connections on E^. This number is a value of the Hilbert function of the toric model of a connected graph 
with no leaves and the first Betti number g. 

By Theorem 14.51 '^6 know that the Hilbert function only depends on the topological invariants of the 
graph. 

Although the model depends on the shape of the underlying trivalent tree, once we restrict ourselves to 
trees with fixed number of leaves, models of all of them are in the same irreducible component of the Hilbert 
scheme of projective varieties with fixed Hilbert polynomial. This was proved by Sturmfels and Xu in [SXlOj . 

Any trivalent graph is made by gluing together tripods, that is graphs with four vertices and three 
edges attached to the central vertex. To construct the toric model we assign to every inner vertex a copy of 
a three-dimensional complex projective space and to every edge we assign an action of the one-dimensional 
complex torus C* on the product of all those P^, which corresponds to gluing two tripods along that edge. 
The model X{0) of the trivalent graph 5 is a geometric invariant theory (GIT) quotient of product of the 
P'^ by the torus defined as a product of the C*'s corresponding to the internal edges. We also translate this 
description into language of projective toric varieties, by writing the model X{Q) as the projective spectrum 
of a semigroup ring C[t{Q)]. The underlying semigroup t(Q) has a clear description in terms of the graph Q. 

Three results of this article generalize our earlier results obtained in |BW07| about binary symmetric 
models of trivalent trees to phylogenetic graph models. First we describe the minimal Z-generators of the 
semigroup t{0) when the graph Q has the first Betti number at most one. We also prove that models of 
graphs with the same discrete invariants are deformation equivalent and lastly that they share the same 
Hilbert function. 
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1.3 Geometric invariant theory. 

We use geometric invariant theory for a normal projective variety X with an action of an algebraic torus T. 
Our main reference is Section 5 and 6 of |BB02| . although the setup we use may seem to be slightly more 
general then the one found in |BB02] . This is because instead of linearizing only with respect to a line bundle 
we allow ample Weil divisors. We will explain the necessary modifications and show how this does not affect 
the basic theory. 

Definition 1.1. A divisor L is an ample Weil divisor if some positive multiple nL is an ample line bundle. 
Given an ample Weil divisor C we have the ring 

oo 

R{X,C) :-0il"(X,O(p£)), 

p=0 

which is the projective coordinate ring of X embedded into a weighted projective space by the linear system 
\C\. This is completely analogous, see jReiOO] . to the standard way of describing embedding of X into a 
projective space in when £ is a very ample line bundle, see |Har77[ Section II. 2]. We discuss those facts 
inOl 

We denote by 

Rp{X,C) := R{X,pC) 

the ring given by a multiple of L. When the divisor C is clear, we write R{X) and Rp{X) an instead of 
R{X,JZ) and R{X,pjZ). 

Lemma 1.2 (Veronese embedding). Let X = Proj R{X, C) be a projective variety with a distinguished ample 
Weil divisor C Then for any positive integer p the inclusion Rp{X, C) ^ R^X, C) induces an isomorphism 

Proj R{X, C) ~ Proj Rp{X, C) 

Proof. See jGro61| Chapter 2, Theorem 2.4.7], and also }Har77| Ex. II 5.13]. □ 

We briefly recall the definitions of quotients simplifying the statements from jBB02| . Let G be a reductive 
algebraic group acting on a variety X . A G-equivariant morphism tt : X ^ Y \s called a semi-geometric 
quotient if 

• the image of every closed orbit is closed, and this property is invariant under base change 

• TT is surjective and images of disjoint, closed orbits are disjoint and this property is invariant under 
base change 

The map tt : X -^Y \s called a good quotient if it is both affine and semi-geometric. A map tt : AT — > y is 
affine if preimage of any afhne subset of Y is affine. 

This definition implies that the good quotient is a categorical quotient, and thus unique |BB02I Rmk.3.1 
and Thm 3.2]. 

To construct the quotient, we need some more definitions. 

Definition 1.3. Let G be a reductive algebraic group acting on X. An ample Weil divisor £ on AT is 
G-linearized if the action of G can be lifted to the projective coordinate ring R{X, C). More precisely, the 
action of G on AT = Proj R{X, C) induced by the lift, is the action of G on X. 
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When £ is a line bundle, a G-linearization is an action oi G on C which is linear on fibers and agrees 
with its action on X. We recall the definition of a >C-semistable point [BB02[ § 6.1], which is independent of 
the multiple of C. 

Definition 1.4. Let C be an ample Weil divisor. A point x € X is £-semistable if there exists a G-invariant 
section of a positive multiple of C that does not vanish at x. We denote by X*" the set of all semistable 
points. The GIT quotient of X by G with respect to the linearization C is 

X//G = X'y/G, 

where X^^ jjG denotes the good quotient of the set of semistable points with respect to £ by the action of 
G. The rational map X X//G is called GIT quotient map and the regular map — > ^l/G is a good 
quotient map. 

Remark 1.5. Both [MFK94| and |BB02| assume that the section in the above definition has an affine 
support, but as we consider only ample C all its sections automatically have affine support. 

The following theorem says that in the above situation a GIT quotient exists and is the projective 
spectrum of the ring of invariants. 

Theorem 1.6. Let G he a reductive group acting on projective varieties X and Y . 

(i) Let C be a G-linearized ample Weil divisor on X. Then 

XII G = Proji?(X, Lf^. 

(ii) Let X ^ Y be a G-equivariant embedding, with Y = Proj i? and X = Proji?//, where I is the 



homogeneous ideal of X in Y. Then I is generated by invariants /i, . . . , /j, the map XHG ^ YUG 
an embedding and the homogeneous ideal of XHG in is also generated by fi, . . . , f. 



IS 
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Proof. In |BB021 §5] the affine quotient is defined for an affine variety X as a spectrum of invariants and 
|BB02[ Thm 5.4] says it is a good quotient. Then in f BB02( §7] there is a characterization of a good quotient 
77 : X — >■ y as a locally affine quotient, that is the map tt is a good quotient if and only if for any open affine 
subset U dY the restriction of the map tt to 7r^^(C/) is an affine quotient. 



We know from jBB021 Thm 6.2.1] that (i) is true for projective space and C — 0(1). Additionally for an 
arbitrary variety X = Proj R and a very ample line bundle C that defines an equivariant embedding into P" 
by |BB02[ §6.3] we have: 

• the restriction of the good quotient morphism tt: (P")*^^ — > (P")**/G is the good quotient morphism 
on X. 

We choose an affine covering lA of the quotient W' HG by sets of the form {/ 0}, where / is G-invariant. 
Then {tt~^{U)\U G U} is an affine G-invariant covering of P", since tt is an affine map. By jBB02[ Thm 5.3] 
and the Hilbert-Nagata theorem (BB02| Thm 5.2] for each U GU we know that (tt^^ {U)nX) HG is a spectrum 
of invariants 

(^-i(t/)nX)//G-Spec(i?[ri]T 
Also each tt^^{U) D X has form {g ~ TT*f\x 7^ 0}. To see that X//G = Proji?*^, we only need to see that 
{^R'^lg^^fj — {R[g^^]^)'-' , which is true since g is G-invariant. Thus (i) holds for X projective and C very 
ample. 
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When C is an ample Weil divisor, we use Lemma [1.21 to replace L by its fc power. More precisely, we 
choose fc such that kC is a very ample line bundle, and since the set of semistable points does not depend on 
the multiple of £ we have 

XIIG = Proj Cf - Proj(i?(X, Cf\ = Proj R{X, Cf 

where the first equality holds, because kL is very ample, the second holds because the action of G preserves 
gradations, and in the third we use Lemma ll.2l 



For the proof of (ii) by Hilbert-Nagata theorem, we know that / is generated by invariants. Since a good 

Y//G is an embedding. The last part of (ii) follows 



quotient is a locally affine quotient the map X//G 
from jBB021 Thm. 5.2]. 



□ 



1.4 Toric varieties in weighted projective spaces. 

The embedding of a projective toric variety X into a projective space is described by a polytope A with 
integral vertices. If we scale the polytope by an integer fc, then we do not change the variety. The resulting 
embedding changes by composing it with fc*'' Veronese embedding. The sum of all positive multiples forms 
a semigroup (or a graded cone). In this situation X — Proj C[fcA]. When the ambient space is a 

weighted projective space the embedding is given by a graded cone with a set of (minimal) Z-generators, 
which are no longer in the first degree. Every section of this cone determined by the grading is a rational 
polytope. 

Definition 1.7. A weighted projective space P(ao, . . . , a„) with weights (ao, . . . , a„) where each G N 
is a positive integer is the GIT quotient of the affine space C"'''^ by the action 

i • (xo, . . . , a;„) = (fo ■ xo, . . . , ■ a;„) 

Again, the non-stable locus is the point G C"+^. 

We can assume that the greatest common divisor of the weights is one — this is by substitution t i— ^ 
^gcd(a(,,...,a„)^ Let us pick integers bo, . . . ,bn such that J27=o bi ■ ai = 1. The weighted projective space has 
the sheaf C'p(l), which corresponds to the ample Weil divisor Y^^^q hi ■ (x^), where the (xi)'s are the divisors 
corresponding to the coordinates. 

Weighted projective spaces are often singular. They have quotient singularities coming from finite abelian 
group actions. 

Definition 1.8. Given a lattice M, we can associate with it in a non-unique way a graded lattice M^"^, 
which equipped with degree map deg : M^^ — 5- Z — the projection to the first coordinate and fits into the 
exact sequence 

^ M ^ MS"- Z ^ 

s 

We also fix a splitting s : Z ^ A/^'' of the exact sequence which is equivalent to a choice of the element in 
the M ~ (1, M) C Mf '. The choice of the splitting s corresponds to a choice of linearization of the action of 
the torus Spec(C[Af^) on itself. 

Definition 1.9. A graded lattice cone r is a rational, convex, polyhedral cone in a graded lattice M^'', 
with all elements having non-negative degree: deg(T) C N, and the zero gradation consists of one element: 
deg~"'^(0) n T = 0. Convex here means that t R_|_ C M^^ M is convex. Equivalently, r is a (saturated) 
sub-semigroup of the free abelian group M^^ with finite set of Z-generators all having positive degrees where 
the neutral element is the only one of degree zero. 
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Proposition 1.10. A toric variety X C P(ao,...,a„) is described by its fan and an ample Weil divisor 
OxO-) or equivalently by an isomorphism X ~ ProjC[r] where t is a graded, rational, convex, polyhedral 
cone in a graded lattice M^^ . Then C[t] = R/I, where R is the homogeneous coordinate ring o/P(ao, . . . ,a„) 
and I is the homogeneous ideal of X . 

Idea of the proof. The correspondence between the ample Weil divisor and the graded cone is the fol- 
lowing. The degree k sections of the sheaf associated to the Weil divisor form the fc*'* section of the cone r. 
To go the other way, we have an isomorphism X ~ ProjC[r] and then the ample Weil divisor is C'x(l) — 
the pull-back of C'p(a(,,...,a„)(l) from the ambient weighted projective space P(ao, . . . , a„). □ 

The choice of the element of M ~ (1, Af) C M^^ in Definition 11.81 is a choice of linearization of the 
action of the torus of X, which extends the action of the torus on itself 

Definition 1.11. A graded product Mf^ Xg Ml^" of the graded lattices Mf" and Ml*" is the fiber product 
over their degree maps, or equivalently the hyperplane degj^ = degj in the product Mf~ x M|'^: 

Mf XgMf >Aff 

dcg2 

Mr ^ z 

Definition 1.12. A graded product ti Xg t2 of graded cones ti and T2 is the fiber product over their 
degree map, or equivalently intersection of the hyperplane degj^ = deg2 with the product cone ri x T2: 

Tl Xg T2 T2 

dcg2 

dcgi ' 

n ^z 

Lemma 1.13. If ti C Mf^ and T2 C M^ are graded rational convex polyhedral cones, and 

Xi=ProjC[ri] X2 = ProjC[r2] 
are corresponding projective toric varieties, then the product of these varieties under the Segre embedding 

Xi X X2 =ProjC[Ti x<,r2] 
corresponds to the graded product of the cones ti and T2 . 
Proof. By definition Xi x X2 under Segre embedding is equal to 



Proj H° (Xi X X2 , Ox, {i)^OxAi))- 

i=0 

For each i we know that the sections of this exterior tensor product are spanned by the product of the z*'' 
graded pieces of the cones 

xX2,OxiWHOx2W) =C - ((nndeg-iW) X (rindeg-i(z))) 

We conclude the lemma by summing the above equality over all z's to get 

00 

^H"{XixX2,OxAi)^OxM^^{^i ^9^2] 

i=Q 

as required. □ 
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1.5 Quotient of a projective toric variety by subtorus. 

We show that the quotient of a projective toric variety X — ProjC[r] by a subtorus is described by an 
appropriate Unear section of r. 

We first recall some facts about toric varieties from |Ful93] . li X ~ Proj C[r], where r C Af is a graded 
cone, then the torus of X is given by an isomorphism T = Hom(M, Z) ®z C*. The lattice N — Hom(Af, Z) 
is the lattice of one-parameter subgroups of T. For any subtorus T' of the torus T there are corresponding 
maps of lattices: the projection M M' of the monomial lattices and the embedding iV' ^ TV of the lattice 
of one-parameter subgroups. 

A linearization of the action of the torus T on X, induces a linearization of the action of any subtorus 
T' C T, by restricting the action. In this setting Theorem 11.61 implies the following result. 

Theorem 1.14. Let t be a graded cone in a lattice M^^ and X = Proj C[r] the corresponding toric variety. 
Let j : T' T 6e a subtorus of the torus T, with : N' ^ N , j* : M M' and id ® j* : MS'' (M')^'' the 
corresponding lattice maps. Then there exists a good quotient and it is equal to 

X//T' = ProjC[Tf' , 

where 

C[t]'^' = C[rn (Z©ker(j* : M Af'))] = C[r n (Z ® 7V'^)]. 
Moreover the quotient X//T' is polarized by Ox//t{^) a natural way. 

The following example shows that we do need ample Weil divisors, not only ample line bundles. 

Example 1.15. Let X be the good GIT-quotient 

TT : p3 X p3 ^ X = p3 X P^)//{C* f 

of the product of two projective three-spaces by an action of three-dimensional torus acting with weights 

"(0 1 1 0) X (0 -1 -1 0)' 
(0 1 -1 0) X (0 0) 
(0 0) X (0 1 -1 0)_ 

linearized with respect to the line bundle C'p3(l) Kl 0^3(1). We will see later that X = X (^) is the model 
of the trivalent graph ? with the first Betti number two with three edges and is a projective toric variety 
by Theorem 11.141 The sheaf Ox(l) = 7r*(C'p3(l) Kl Op3(l)) is not a locally free 0x -module because the 
associated divisor is not Cartier. To verify it we can use a computer algebra system, for example magma 
|BCP97) as follows. Since any divisor on a toric variety is linearly equivalent to a T-invariant divisor, we 
identify a divisor with a corresponding Z-combination of primitive elements of the rays of the fan. Thus we 
only need to check if the Z-combination corresponding to Ox(l) yields a piecewise linear function on the fan, 
which by |Oda88j is equivalent to our T-invariant Weil divisor being Cartier. 

2 Phylogenetic models on trivalent graphs 
2.1 Trivalent graphs. 

We define topological invariants of trivalent graphs and show any two graphs with the same invariants are 
equivalent by applying appropriate mutations, which we introduced in [BW07) . We do not assume that our 
graphs are connected. 
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Definition 2.1. A graph Q is set V of vertices and set £ of edges together with the unordered boundary 
map :£ v®2^ where V®^ is the set of unordered pairs of vertices. We write (e) = {9i(e), 92(e)} and say 
that V is an end of the edge e if w S (e). A vertex incident to exactly one edge is a leaf. The set of leaves 
is denoted by C and the number of leaves by n. If a vertex is not a leaf, it is called an inner vertex. An 
edge incident to a leaf is a petiole and V is the set of petioles. We write comp Q for the set of connected 
components of the graph and | compCyj for the number of components. We denote by g for the first Betti 
numbers of graph, which is the rank of the first group homology of the graph viewed as a CW-complex. A 
graph is trivalent if every inner vertex has valency three. Valency of a vertex ti is the number of connected 
components of a sufficiently small neighborhood of v with v removed. A trivalent graph with no cycles is a 
trivalent tree. 

When discussing more then one graph instead of V, £, n,... we will write V{G), £{G), n{Q), etc. 

Remark 2.2. Our graphs are not oriented, nevertheless we write 9i(e) and 92(e) for the vertices adjacent 
to the edge e. This makes it easier to talk about "the other end of e" . 

We call the unique trivalent tree with a single inner vertex and three leaves the elementary tripod. 
It has three edges ei, 62, es. Any trivalent graph is built of elementary tripods in the following way: given a 
trivalent graph Q and any inner vertex v S V(^) we pick a copy of elementary tripod J^y ~ and an map 
ill : J^v — ^ G which sends the central vertex of J^^ to v and locally near i~^(v) is an embedding. We present 
the graph Q as a, disjoint union of the graphs J^^ with appropriate identification of edges 

U -^'^/{*ai^(e)(e) «a2V)(e)}ee£\7' (2.3) 
vev\c 

This construction mirrors how the model of the graph is constructed, as we will see in Definition 12.121 

Example 2.4. On the Figure[T]we give an example of the above presentation of a trivalent graph for a graph 
with the first Betti number one and two leaves. 




Figure 1: Building a trivalent graph from tripods 



Lemma 2.5. In any trivalent graph with n leaves and first Betti number g the following holds 
ft) |V|,|^| >n, 

(ii) 2\£\ = 3|V| - 2n, 

(Hi) \V\ - \£\ = \compQ\ — g. 
Thus, any three of the numbers \V\, \£\, n, g, \ coTapQ\ determine the other two. 



Proof. To prove (ii) let us count pairs of consisting of a vertex and an adjacent edge. On one hand we 
will count every edge twice. On the other hand every inner vertex has three incident edges so we have 
3(|V| — n) pairs and another n pairs come from leaves which totals to 3|V| — 2n. Equation (iii) counts the 
Euler characteristic | compt?| — g of the graph. □ 
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Here we introduce operations of gluing two leaves of a graph, cutting an edge into two new edges and 
taking a disjoint sum of two graphs. 

Definition 2.6. We will use the following three constructions of trivalent graphs. 

• Q \JQ' is the disjoint sum of the given graphs Q and Q' . 

• GI^D is the graph obtained from a given graph Q with two distinguished leaves li,l2 £ ^{G) by gluing 
the two leaves li and I2, or more precisely by removing li and I2 and identifying the edge incident 
to li with the edge incident to l2- 

• Gi-^^Q'i a graft of given graphs G and G' each with a distinguished leaf. Figure [2] is a schematic picture 
of this construction. 



• G"^ is the graph obtained from the given graph G by cutting an internal edge e G ^{G) \ ^{G)- More 
precisely we replace e by two new edges ei and 62 with 9i(ei) := 9i(e) and 81(62) := 92(e). There are 
two new leaves in G'^, which are the free ends of the new edges 82(61) and 82(62)- 

Definition 2.7. An edge e G f is called a cycle edge if it is not a petiole and removing it does not disconnect 
the connected component of the graph that contains e. An edge e G £ is called cycle leg if it is incident to 
a cycle edge but it is not a cycle edge. A vertex u G V is called cycle vertex if it is an end of a cycle edge. 
We draw example of those on Figure [3] 




Figure 2: Graft of two graphs 



The new graph can be written as 



Gi*G[, = {{(GuG')\p)uG')[p 



cycle legs 




cycle edge 



cycle vertex 



Figure 3: Cycle edge, cycle leg and cycle vertex of a graph 
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A path is a sequence of distinct edges eo, . . . , with d2{&i) — 9i(ei+i) for all i G {0, . . . , m — 1}, moreover 
c?i(eo) and d2{em) are both leaves or they are either equal. In the latter case, the path is called a cycle. 
Paths are disjoint if they have no common vertices. A network is a union of disjoint paths. For consistency 
we say that the empty set is also a network. A cycle is a minimal sequence of cycle edges. A cycle of length 
one is a loop. On Figure |4] we draw examples for each those sequences. 





Figure 4: A path, a cycle and a network containing a loop 



A graph Q is called a polygon graph if it has 2k edges of which k form the only cycle of Q and the remaining 
k edges are cycle legs. If Q is any graph, e G V{Q) a non-cycle edge and after cutting e we get a decomposition 
Q"^ — QqU Qi where Qi is a tree then we call Qi a pendant tree. 

Figure [5] shows the three trivalent trees with one internal edge e and four labeled leaves. 




Figure 5: The three trees with four labeled leaves 



Let Q he a trivalent graph and let e be an internal edge which is not a loop (the ends of e are not 
identified). Then a neighborhood of e in is a trivalent tree with four leaves. 

Definition 2.8. A mutation of a trivalent graph Q along an edge e is a graph Q' which is obtained from 
Q by removing a neighborhood of e (which is a trivalent tree with four leaves) and replacing it by one of the 
other trivalent trees from Figure [5] Two graphs are mutation-equivalent if they may be transformed by 
sequence of mutations into each other. 

Remark 2.9. By definition it is not possible to mutate along an edge that forms a cycle of length one. 
However, for longer cycles, mutations are possible, and one shortens the length of the cycle. Figure [6] shows 
this fenomena for Hammock graph mutation equivalent to LittleMan. 




Figure 6: Mutation along cycle edge shortens a cycle 
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Lemma 2.10. Suppose edges {ei,...,ek} form a cycle in the graph Q and assume k > 1. Then, for any 
i € {1, . . . , fc} one of the two mutations along shortens the cycle {ei, . . . , 6^} by one in the resulting graph, 
i.e. {ei, . . . , Ci-i, Bi+i, . . . , Cfc} is a cycle in the new graph. 

A caterpillar is a trivalent tree, which after removing all leaves and petioles becomes a string of edges 
as shown on Figure [T] 




Figure 7: Caterpillar tree 



Lemma 2.11. Let Qi and Q2 he connected, trivalent graphs both with n leaves and first Betti numbers g. 
Then they are mutation-equivalent. Moreover, for any ordered subsets Si of cycle edges of Qi and S2 of cycle 
edges of Q2, of the same size, both with the property that removing Si from Qi does not disconnect the graph, 
we can find a sequence of mutations that avoid the edges from Si (S2) and sends i*^ edge of Si to the i^^ 
edge 0/6*2. Also, any mutation sends a leaf of Qi to a leaf of Q2. 

Proof. Let Q be connected, trivalent graph with n leaves and the first Betti number g and S a subset of 
cycle edges as above. We will prove that Q is mutation-equivalent to a trivalent graph obtained by attaching 
g cycles of length one to a caterpillar tree with n + g leaves. We will choose mutations so that they will 
satisfy the required property. 

Step 1. We proceed by induction on the set S and the first Betti number of Q. For an edge e G we can 
find a cycle that contains no other elements of S. This is because after removing all edges from S the graph 
Q is connected, so there is a path 7 from 9i(e) to 92(e), which together with e form the required cycle. By 
repeatedly using Lemma l2.10l we reduce the length of this cycle to one, by performing mutations along edges 
from 7. In the new graph the edge e forms a loop. We can consider this graph with e removed, reducing 
both the size of S and g. If \S\ < g, then we repeat the above g — \S\ times starting from any cycle edge, 
which is not a loop. 

After repeating this procedure g times, we get a tree with g loops (all edges from S are among them) 
attached to some leaves. We can assume that this tree is a caterpillar, as we know from [BWO?) Lem. 2.18], 
that any trivalent tree is mutation-equivalent to a caterpillar with the same number of leaves. 

Step 2. We observe that it does not matter to which leaves the cycles are attached, we can move a cycle 
from a leaf to any another leaf. 

In Figure [8] we illustrate those two steps. 




Figure 8: Every graph is mutation-equivalent to a caterpillar graph 
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The last claim follows simply form the definition: mutation maps an inner edge to an inner edge, and a leaf 
to a leaf. □ 



2.2 Definition as quotient. 

Given a not necessarily connected trivalent graph Q, we construct a toric variety X{Q), generalizing the 
binary symmetric model of trivalent tree from B W07i . 

As we have already explained, see equation (j2.3p . any trivalent graph is the union of |V| — n elementary 
tripods with some edges identified. To define the variety X{Q) we replace each elementary tripod J^^ with 
P^, union with product, and the edge identification with a quotient by an action of a one-parameter torus. 

Definition 2.12. Let Q he a, trivalent graph. To an inner vertex v & V\C we associate projective space 
with coordinates X0,xl2,Xi^,X23- To any edge e G £" we associate an action of C* on with weights 
and 1 as follows: 



t ■ x"g if the index of ^[e) G {ei, 62, 63} belongs to the set 5", 



x^g otherwise. 



So we have an action of a three-dimensional torus on P = Proj <C[xq|,x\2^ 0:23] with weights 

110' 
10 1 
11 

In other words if for example it, (62) = e then acts with weight 1 on x\2 and 2:23 and with weight on x^ 
and x\^. This action extends to an action A^ on Iluevv^^u which is non-trivial only if v is an end of the 
edge e. Thus, for any internal edge e G £ \ P, we can define a C*-action A^^j^-j x "-^^^(e) on 

n (2-13) 

vt^V\C 

to be the product action of the action A^^ on P|^ and the action A^^ ^^^j with opposite weights on P|^ . 
We define the phylogenetic model of a trivalent graph Q to be the good quotient: 

X{Q):={ n / n (Aa.(e)X-A^,(e)) (2.14) 

v&V\C e££\V 

of the toric variety n^evvc'^? ^ subtorus of dimension \£\ — \C\ of the torus. The subtorus by which we 
are dividing is a product of all the C* 's over all internal edges of the graph Q and the linearized line bundle 
is My^vi^g)Opa {l). By Theorem II. 141 X{Q) is toric as it is the quotient of a toric variety by a subtorus. 

Remark 2.15. In Definition (|2.14p the choice that we made defining the action of the torus (C*)^^'^ only 
depends on the choice of coordinates of the torus. If we choose different orientation of the edge e, then the 
two C*-actions -^^^(e) x ^-^^^(e) and —-^^^(e) x Ag^^g-| differ by composing with f i— > i. 

Remark 2.16. Let I be a leaf of a graph Q and e the adjacent petiole. The action Xf descends to a non-trivial 
action on the quotient variety X{Q) and is denoted by A'. For a subset 5 C £ of the leaves of Q of cardinality 
k we have an action of a A:-dimensional torus T{S) — a product of the corresponding A''s. 
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If we set k = \V\ — n the number of inner vertices, we can rewrite |V| — \£\ — |compt/| — g using 
Lemnla l2.^^ii)| to get k = (2g — 2| compel) + n. Observe that 

• A: is the number of the P^'s in (|2?T3l) 

• g is the first Bctti number. 

On the other hand, \£\ — n = k + g — \ compel is the number of inner edges which is the dimension of the 
torus that we divide by in (I2.14p . We get a variety of dimension 

dimX{g) = 3fc - (1^1 - n) = 3g - 3| compel +2n = \£\. 

2.3 Lattice and cone. 

Given a trivalent graph Q we construct the toric data that allows to recover its toric model. The graded 
lattice, denoted by Af^'', as well as the graded cone in it have both rank one bigger than the dimension of 
the model X{Q), and the latter is equal to the number of edges £{9)- 

Definition 2.17. Given a graph Q let Z£ = Z • e be the lattice spanned by £, and ^ = Hom(Zf , Z) be 
its dual. Elements of the lattice Zf are formal linear combinations of the edges, forming the standard basis 
of Zf. The dual lattice Zf^ comes with the dual basis {e*}e££- We identify vertices of the graph Q with 
certain elements of Zf ^ : 

v = Y,e*. (2.18) 

We also define AI = {u e Z£ -.Vv e V v{u) e 2Z} and its dual N = Hom(M, Z). Then the graded lattice of 
the graph is 

MS"- = Z © M, 

with the degree map 

deg : MS"" = Z © M ^ Z, 
which is the projection to the first summand. The degree of w e M^" is deg(u'). 

If there is more then one graph in question we will write M{Q) for M^^{Q) etc. 

Let us use the following notation for the elements of the lattice Zf ^ dual to the edges meeting at the 
vertex v 

ay := («t,(ei)) , 6„ := {iv{e2)) , Cy {ivies)) , 

where {ei, 62, 63} are the edges of and : ^ Q is, as before, a map which is locally an embedding and 
sends the central vertex of the to v — an inner vertex of G- Whenever we use this notation we have a 
fixed presentation as in (|2.3p . 

Given an element uj in either Z£, M or Af^*", each of ay,bv,Cy G Zf ^ measures the coefficient of w at an 
edge incident to v. Then (|2.18l) becomes 

Definition 2.19. The degree of a; G A/^'' at a vertex v G V{Q) is 

deg„(w) := i • {ay{Lu) + by{uj) + Cy{u))). 

The minimal degree of cj is 

degmi„(w) := max{degjw)}, 

i>ev 

where ttm : M is the projection to the second summand. 



14 



The name minimal degree will be clear after we define the cone t(Q). 

We identify paths and networks in Q with elements of the lattices M and M^'^, by replacing union with 
sum in the group 

Definition 2.20. A network in the graded lattice M^^ is a pair uj = (l,a) e M^"" where a e M is a 
network. 

Lemma 2.21. An element of the lattice M is represented by a labeling of the edges of Q with integers so that 
the sum at any vertex is even. Thus the lattice M C is generated by 

(i) networks 

(ii) {2e I e e £}. 

Proof. Let w G M. By using generators of the second type, we can assume that < a^{uj), by{uj), c„(w) < 1 
for any vertex v. This implies that a^(cj)+6„(a;)+c„(w) equals or 2 and either two among a^(w), 6^(cj), Cv{uj) 
are one or all are zero. This implies that cj is a network, since it corresponds to a disjoint union of path. A 
path goes through a vertex v means in terms of ai,(cij), 6„(a;), Cv{io) that exactly two of them are one. □ 

We define the cone t(Q) of the graph, which is the semigroup defining the model of the graph as projective 
spectrum of the semigroup ring, as we will see in Theorem 12.241 

Definition 2.22. For a graph Q we define its cone r = t{Q) C M^^ as the set of cj G M^^ which satisfy 
following inequalities: 

(i) ay{iAj),by{uj),Cy{ul) >0, 

(ii) for any vertex w e V triangle inequalities hold 

|at,(a;) — by{u!)\ < Cy{Lu) < ay{uj) + bv{uj), and 

(iii) deg(a;) > deg,„i„(a;). 

Remark 2.23. To explain the name minimal degree degjjj;„(a;), note that for any lo in the cone we have the 
following equality 

deg™„(w) = min{deg(w') : tTj^{uj') = TT^iuj)}. 



Proof of the remark. By part (iii) of the Definition 12.221 of t for any w' G t satisfying TT-g-{u!') = n^{uj) 
we have 

deg(a;') > deg^i„(w') = deg,„i„(a;) 
since by definition of degj„jj-|(aj) only depends on ■n^{uj). This means 

deg„i„(a;) < min{deg(tj') : 7rjj(a;') = 7rjj(a;)}. 

To prove the equality we will find w' e r with deg(a;') = degjjji,j(a;). Let us write uj = (deg(a;),a) as it is an 
element of M^^ . Recah that deg„i„(a;) is the maximum of deg„(aj) = 5 • ( ){w) over ah vertices v 

of the graph. Thus uj' = {deg^^^{u}) , a) is in the cone r and has the required degree in M^*". □ 

Theorem 2.24. The variety X{Q) is isomorphic to the toric variety Proj C[T((y)]. 
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Proof. To see this we first observe that each in Definition 12.121 of the model X{Q) can be written as 
Proj C[t(/L,,u)], where each cone 

t{Xv) = conv(0000, 1000, 1110, 1101, 1011) 

is a cone over a tetrahedron and is clearly defined by the required inequalities. Next, taking the product of 
corresponds by Lemma [1.131 to taking a graded product of cones. Thus the product cone is defined by 
required inequalities. Lastly we use the description of the quotient of a toric variety by a subtorus of the 
torus in Theorem ll.141 The subtorus in Definition 12. 121 is a product of the C* over all internal edges. Taking 
the quotient with respect to such a torus corresponds to cutting the cone with the hyperplane of the type 
ddiie) = bo2{e), which preserves the inequalities. □ 

2.4 Z-generators of the cone r{Q) 

Knowing that the model X{Q) is the projective spectrum of the semigroup algebra of t{Q) means that it 
is a subvariety of a weighted projective space with weights equal to the degrees of the chosen generators. 
When Q = T is a tree the cone is generated in degree 1 so the embedding is into a (straight) projective space 
pfc _ p^j^^ . . . , 1). In this case, by |BW07] . we already know all about this cone, see Proposition 12.321 below. 
It is represented by its degree 1 section — a normal lattice polytope A(T) C 1 x M, whose vertices span the 
cone t(T) and the lattice points generate the semigroup. Our goal is to show that for graphs with the first 
Betti number one, the semigroup is generated in degrees 1 and 2. 

In order to describe Z-generators of the cone t{Q) we will express elements of t{Q) in terms of Q. We 
also decompose the graph Q into smaller graphs for which Z-generators of the corresponding cones are easier 
to find. 

We explain that any element uj of the cone t{Q) locally decomposes into paths. In the graph there 
are three non-empty paths, each consisting of two edges. Let us denote them by: 

x:^e2 + e3, y:=ei + e3, z := ei + 62, 

where {ei, 62, 63} are edges of J^. 

For an arbitrary Q, we know that given an element uj G t(Q) of the cone and a vertex v G V{Q), the 
numbers ay{uj),by(u)),Cy(LL)) satisfy the triangle inequalities and their sum is even. This allows us to, locally 
at V, rewrite w as sum of paths x, y and z. The picture of this decomposition is drawn in Figure |9l 




Figure 9: Local paths around a vertex 



Our aim is to find the Z-generators by understanding how the graph Q was built from smaller pieces. Each 
of the operations in Definition 12.61 has a corresponding operation on lattices and cones. By Definition 12.121 
the model of a disjoint sum of graphs is the product of the models, so the underlying cone is the graded 
product of corresponding cones. 



16 



Lemma 2.25. Let Qi and Q2 be two trivalent graphs then 

Ms^iGi U 02) = M9'-{gi) Xg M9^{g2), and 

T{giug2) = T{gi) XgT{g2). 

In the definition of X(g) we talce a quotient by a torus corresponding to tlie set of inner edges. In 
other words we have translated the operation of gluing of two leaves of a graph g into taking a quotient by 
appropriate C*-action of the model X{g). The following observation says that if we choose to glue some 
pairs of leaves first and then the rest of the pairs it does not matter how we partition the set of pairs of leaves 
or which order we choose. In all cases the resulting variety is the same. 

Lemma 2.26. Let two reductive, commutative groups Hi and H2 and their direct sum Hi ® H2 act on a 
projective variety X. Suppose all those actions are linearized with respect to some ample Weil divisor C 
Then 

X//{Hi © H2) = {X//Hi)//H2 = {X//H2)//Hi, 

where the semistable points on X are taken with respect to C and on quotients of X with respect to the 
push-forward of C. 

We have identified the vertices of g with elements of ^, see ()2.18|) . We observed in Remark 12.161 that 
a leaf I yields a C*-action A' on X{g). Given two leaves li and I2 of g, by definition we have 

x{gliz))^x{g)//{x'^ x-x'^). 

In terms of toric geometry this quotient corresponds to the intersection of t(C?) with the kernel of li — I2, 
where we treat li and I2 as elements of the lattice (M^'')^. Thus the following lemma is a consequence of 
Theorem 01 

Lemma 2.27. Let li and I2 be distinct leaves of g. Then 

M9^{g\p) = M9-{g) n ker(/i-Z2) 

T(5|p) - T{g) n keiih-h). 

The reverse operation on graphs is to cut an edge into two new edges. In the next lemma we explain 
how this is refiected on the cones and lattices. 

Lemma 2.28. Let g be a trivalent graph and e £ \ 'P an internal edge. As before g'^ is the graph obtained 
from g by cutting the edge e. Then there are natural embeddings of the cones and graded lattices: 

: MS'-(^) ^ M3-{g^) : t(^) ^ T{g') 

Proof. Let 

Zr(g)= Z-e' 

e'e£{g)\{e} 

be the lattice spanned by all other edges. We can decompose the lattices {Z£)9^{g) and {ZE)^^ {g"^) 

Z£{g) =Z(BZ£''{g)(BZ-e, 
ZSig") = Z ® ZS^ig) © Z ■ ei © Z • ea. 
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Thus we can embed the lattices identifying the first two summands and taking a diagonal embedding of the 
third one: := idz®idz£= ©A^ : {'Z£)^^{0) ^ {Z£)^^{0'^). To complete the proof we check that i*^ restricts 
to the lattices with parity condition and to the cones. 



A A 



T{g^ 



To see it we only need to check the parity condition about the ends (3i(e) and 92(e) of the edge e that 
we cut. By definition e*(i^(e)) = e2{i^{e)), so for any uj e M^"^ we have di{e){uj) = 5i(ei)(i'^(w)) and 
d2{e){uj) = 92(e2)(«'^(i^)). In the same way i'^ preserves all the inequalities defining cones T{g) and T{g'^). □ 

When the edge e is not a cycle edge, the graph g'^ is not connected. We write = t/i U C/2, where 5i 
(respectively ty2) is the part containing ei (respectively 62). Then we have a projection ttJ (respectively 
of lattices 

ttI : M^-ig") = MS^igi) X M3-ig2) ^ AF"(ei) 

which restricts to a projection of cones. For a non-cycle edge e G £ we denote by pi (respectively P2) the 
composition pi ^ o i'^ of the above defined maps. 

Remark 2.29. When e € £ is not a cycle edge we write g^ = giU 02- Then the cone t(^) is the following 
fiber product of the cones Pi(t(^)) = T((yi) and P2{T{g)) = T{g2)- The same is true for the lattice M9^{g). 



r{Q)- 

P2 



r{Qi) 

dcg ©ej 



^ dog ffieS 

T{g2) — 



Ms^{g) — ^M9''(ai) 
P2 

M9'-{g2) i>z©z 



Now we turn our attention to our main task of finding the Z-generators of the cone T{g). 
Lemma 2.30. For any graph g the set of degree 1 integer points of cone T{g) is equal to the set of networks. 
Proof. If G T is a point in the cone of degree 1, then for any vertex u e V, 

1 = deg(w) > deg„,i„(cj) = max{deg„(a;)} > deg„(a;) > 0. 

By definition deg^ = Xy +yv + > so exactly one of Xv{uj), y«(a;), Zv(u!) equals one and the other two are 
zero or all are zero. Equivalently exactly two of 0^,(0;), by(uj), Cv{uj) are one, and the third one is zero, or all 
are zero. This means that w is a network. □ 

Corollary 2.31. All networks are among the minimal "L-generators of the cone T{g). 

In fact when the graph in question is a tree these are the Z-generators. 



Proposition 2.32 ([BW07, §2.1]). If g is a trivalent tree, then T{g) is generated in degree 1. Moreover, the 
generators are exactly networks of paths, which in this case are determined by their values on the leaves. Thus 
a generator of t(C/) is identified with a sequence of 's and 1 's of length n{g) with even number of 1 's. As a 
consequence a model of a trivalent tree with n leaves comes with an embedding into projective space P'^ 
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Let Q he a, graph with the first Betti number one. We will describe the generators of the semigroup t{Q) 
in this case. We cut all the cycle legs li, . . . Jk of Q and write C/'i> - '''= = Qq U Qi U . . . \J Qk, where Qo is a 
polygon graph and thus Gi, ■ ■ ■ ,Gk are pendant trees (see Definition 12. 7|) . Thus any element uj € t{Q) has a 
hft UJ G t{Qq U . . . U Qk) = t{Go) Xg t{Qi) Xg . . . Xg T{Qk) and components w = (wq, wi, . . . , tffe), which can 
be written LxJi — 

Theorem 2.33. Let Q be a trivalent graph with the first Betti number exactly one. Any element lo G M^''{Q) 
is a minimal Z-generator of t{Q) if and only if it satisfies one of the following conditions 

(i) Lo has degree 1 and uj is a network, or 

(ii) ijj has degree 2, and satisfies the following three conditions determining Wq 

e*{uj) — I, for all cycle edges e G £ \ £ 
e*{uj) = 2, for an odd number of cycle legs, 
e*{uj) = 0, for the remaining cycle legs. 
Each of the remaining components LOi G T(Qi) is an element of degree at most two. 

We postpone the proof until we prepare for it with some lemmas. The idea of the proof is to use 
Lemma 12.281 in order to remove the pendant trees Qi, . ■ ■ ,Qk and work only with the polygon graph ^o- 
Lemma [2.371 describes all the degree 2 points of the cone of a polygon graph and distinguishes the generators 
among them. 

Example 2.34. In Figure [10] we illustrate the generators of the cone associated to the graph LittleMan (one 
of the two graphs with one cycle and two leaves). The first four are of degree 1, the remaining three are of 
degree 2. 




Figure 10: Generators of the cone of the LittleMan graph 



Lemma 2.35 (Decomposition propagates to pendant trees). Let Q be any trivalent graph and w G t any 
cone element. Let us also fix a non-cycle edge e £ \ V such that Q'^ is a tree, where Q'^ = Gi ^ G2- Then 
any decomposition of p\{lj) lifts to a decomposition ofuj. 

Proof. First note that both pf and preserve the degree, so an element w of degree d in t{Q) yields 
wi S t{Qi) and uj2 G '''(&) both of degree d. The semigroup of a tree is generated by networks, which are 
degree 1 elements, see Theorem l2.32l This means that p|(a;) is a sum of degree 1 elements. Thus if pf (w) can 
be decomposed, then the same decomposition works for uj by choosing appropriate grouping of the summands 
of p|(w), because the degrees are preserved. □ 

Corollary 2.36. In the proof of Theorem \2.SS\ we can assume that the graph is a polygon graph. 
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Proof. A graph with one cycle is a polygon graph with a tree attached to each cycle leg li, . . . ,lk- We cut 
all the cycle legs to obtain k + 1 pieces of the graph Q: a polygon graph Qq and k trees: Qi,. . . , Qk- We 
denote by po the composition of projections for each leg that we cut po — pjj o ■ ■ ■ o p^^ . Iteratively using 
Lemma l2.35l to decompose po{u}) we decompose uj. □ 

Lemma 2.37 (Degree 2 elements of the cone). Let Q be a graph with exactly one cycle. Any degree 2 element 



of t{Q) except those in Theorem \2.3S\ (ii) is a sum of two networks 



Proof. By Corollary 12.361 we may assume Q is a. polygon graph. Let uj G t{Q) a degree 2 element. The 
coefficient e*{uj) of u) on any edge e G f is 0, 1 or 2. We denote by 

L0e\-P ■.^{e*{Lo)\e^£\V} 

the set of coefficients of w on the cycle edges. We distinguish between four types of uj based on ujgyp. For all 
but one we decompose cj as a sum of two networks. 

If G tos\Vy there exists a cycle edge e G £{Q) with e*{uj) = 0. We can cut it with no harm to uj, since 
i'^{uj) G t{Q'^) is a degree 2 element in a cone of the trivalent tree G'^, so it can be decomposed into a sum 
of degree 1 elements. This decomposition can be lifted to t{Q), as we assumed e*{uj) — 0. On Figure [TT] we 
show an example of this situation. 




Figure 11: Decomposition of a degree two element with a weight zero cycle edge 



The second case is when uj£\-p = {2}, that is e*(a;) is 2 on all cycle edges. As oj has degree 2, we know that 
deg^,(w) < 2 and as a consequence: 




2 if e is a cycle edge 

otherwise, i.e. e is a cycle leg. 



Thus Ld is twice the network consisting of all the cycle edges, as on the example on Figure [T2l 




Figure 12: Decomposition of a degree two element with all cycle edges of weight two 
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For the last two cases we let Zi, . . . , Ip be the set of all cycle legs with ^ ordered anticlockwise 

and indexed by elements of Zp. Also we denote hy i -r^ j the path that starts at k and goes anticlockwise 
along the intermediate cycle edges to Ij and ends there. 

In the third case uJe\'p = {1} and we will show that lo can be decomposed into a sum of two networks 
if and only if the number of cycle legs for which uj has a coefficient 2 is even. Suppose we have decomposed 
= wi + W2 into a sum of networks and neither wi nor 0^2 contains the path consisting of all cycle edges 
(in which case the other one would be an empty network). Then both uji and UJ2 contain a positive even 
number of legs. Any path in wi (respectively UJ2) is of the type i -r\ i + 1, from the cycle leg k to the 
next one. If the end was not the next one, there would be a path in i02 (respectively cji) containing an 
intermediate leg and as a result there would be a cycle edge e with w(e) = 2, but w^^-p — {1}. Moreover, 
both (Ml and UJ2 contain all non-zero cycle legs, each with value 1, because they are both of degree 1 and 
in the cone we have deg(a;) > deg^j„(w). This is only possible when the number p of non-zero cycle legs 
is even, and in that case we have the obvious decomposition with uJi = h ^ h + ■ • ■ + ip-i ^ ip and 
0^2 = *2 ^ *3 + • • ■ -I- jp ^ «i . Otherwise is a generator. Examples of both these situations are drawn 
on Figure [T51 




Figure 13: Typical generator of degree two and a decomposition of a degree two element with 
all cycle edges of weight one 

In the last case w^y-p — {1, 2}. When l*{uj) = 2 we call / a two-leg. Denote hy k^,. . . , the subsequence of 
two-legs, numbered in such way that traveling along the cycle anticlockwise from to there is a one-leg, 
provided that uj has some one-legs. We observe that the number of one-legs between two consecutive two-legs 
is always even. This is best explained by drawing the picture from Figure 1141 




two— leg one — leg next two — leg 



Figure 14: Element of degree two has even number of one-legs between consecutive two-legs 



where the arcs are our z^'s introduced in[9l To produce an element of the lattice, the two on the same 

edge need to share the same value. If there would be only zero-legs where the first dots are, the local paths 
would not agree on some cycle edge. 

The decomposition w — (1, uji) + (1, UJ2) depends on the parity of g, which is the number of two-legs. We 
first work in the case g = 2r is even, the odd case uses the same idea with small modifications. The Figure [15] 
visualizes how the decomposition is constructed in the case q = 4. 
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Figure 15: Decomposition of a degree two element with even number of two-legs 



First we place all paths between two consecutive two-legs starting at an leg with an even (respectively odd) 
index into wi (respectively W2). Then, to take care of the one-legs, we add paths between consecutive one-legs 
lying between U^. and ^iaj+i foi' some j S {1, . . . ,p} to oji. Thus we get 

n^i2 + (12 + 1)^(12 + 2) +...+ (43-2)^(43-1)+ 

i3-r\i4 + (14 + 1)^(14 + 2) +...+ (45-2)^(45-1) + 

Wl = . 

«2r-l ^ 42r + (42r + 1) ^ (42r + 2) +...+ (il — 2) -/^ (il - 1) 

Similarly we add paths between consecutive one-legs lying between Ziaj-i ^^'^ ^i2j for some j G {1, . . . ,p} to 
UJ2. So we can write 

i2 ^ 43 + (43 + 1) -r^ (43 + 2) + . . . + (44 - 2) ^ (44 - 1)+ 

44^15 + (45 + 1)^(45+2) +...+ (46-2)^(46-1) + 
42r^4l + (41 + 1)^(41+2) +...+ (42-2)^(42-1) 

Clearly those paths in uji (resp. W2) are disjoint, so both are networks and by construction they yield a 
decomposition of our w. 

When the number q of two-legs is odd we need to adjust the above decomposition. Again, we draw an 
example on Figure [16] for q — 3 




Figure 16: Decomposition of a degree two element with odd number of two-legs 



As there is an odd number of two-legs, we will use two consecutive one-legs lying between li^ and to make 
up for the missing two-leg, and we proceed as before to get 

41 ^ 42 + (42 + 1) ^ (42 + 2) + . . . + (43 - 2) ^ (43 - 1) + 

43 ^ 44 + (44 + 1) ^ (44 + 2) + . . . + (45 - 2) ^ (45 - 1)+ 

uji = : 

42r-l ■r\ 42r + (42r + 1) ^ (42r + 2) + . . . + (42r + l — 2) -r\ (42r + l — 1) 

42r + l ^41 — 1 
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and 

42 ^ 43 + (is + 1) ^ (is + 2) + . . . + {ii - 2) ^ {ii - 1)+ 

14 ^ is + (is + 1) -/^ (is + 2) + . . . + (ig — 2) -/^ (ig - 1)+ 

LJ2 = : 

i2r ■r\ i2r + l + («2r + l + 1) ^ (i2r+l + 2) + . . . + (ii - 4) ^ (il — 3) 

ii — 2 ^ ii 

This ends the proof of the lemma about the decomposable degree two elements of the cone. □ 

Proof of Theorem 12.331 The proof yields an algorithm for decomposing an arbitrary element ui of the 
cone t{G) into a sum of degree 1 and 2 generators. By Corollary 12.361 we only need to prove the theorem 
when ^ is a polygon graph. First we fix an orientation of the cycle of Q and we call it anticlockwise in order 
to think of a planar embedding of the graph. Let lu G t{Q). We will find an element /i of degree at most 2, 
such that u ~ jj, G t{Q). Let w be a vertex and the cycle leg attached it. We choose an map iv ■ J-^ ^ G so 
that the edge 63 is mapped to ly, the edge 62 is mapped to the edge which points anticlockwise from v, and 
the edge ei is mapped to the edge which points clockwise from v. In this notation the coefficient of uj at ly 
is measured by Cy and coefficient of the edge anticlockwise (respectively clockwise) from v is measured by ay 
(respectively 61,). We will also use the local paths Xv,yv,Zy defined in ([9]). 
Now we are ready to describe the algorithm to find /i. 

Step 1. If there is a cycle edge e with e*{uj) — 0, we cut e and obtain the graph which is a trivalent tree. 
Thus p'^{ijj) is a sum of networks of path and this decomposition lifts to graph G- 

Step 2. Otherwise e*(w) > 1 on all cycle edges e. We set ^ to have value 1 on every cycle edge. Equivalently 
is defined by setting at every cycle vertex Zy{ii) — 1, Xy{^) = 0, yy{ii) = 0. Now if a; — /j, £ t{G) we are 
done. Otherwise w — ^ fails one of the inequalities defining t{G)- It is not the one with degrees, since for 
each vertex z; G ^ we have deg^,(/i) — 1, which implies deg^^^{u] — /i) < deg(a;) — deg(/j,) = deg(ti; — /x). 
Clearly all coefficients of w — are positive. Thus at some vertex w G V our uj — fi fails one of the 
triangle inequalities. 

Step 3. We will adjust fi to fix the triangle inequalities for uj — 11. If a triangle inequality for uj ~ fi a,t v fails, 
then this is because Zy^u) — 0. In such a case we set ^{ly) — 2, which will not make any coefficient 
of w — /i negative provided Cy^cd) > 2. But since ay{Ltj),by{uj) > 1 and Zyito) = we must have 
Xy{u!),yy(Lo) > I. This implies Cy{u!) > 2 as required. In terms of Xy,yy,Zy we have decreased Zy{iJ,) 
by one and increased both x^(/i) and ?/«(/i) by one. 

Step 4. We need a last adjustment on /i to assure the additivity of degree where it is attained, i.e. for any v 
such that deg„(a;) — deg(w) we need deg^(^) = 2 since deg(/i) — 2. This is to ensure uj — ^ E t. We 
call V degree deficient vertex if deg^(a;) = deg(w) and deg„(/i) = 1. 

If V is degree deficient and in addition Xy{uj) > and yy{uj) > we set Xy{p) = yy{fJ,) = 1 and 

Zy{^l) = 0. 

If V is degree deficient and both Xy{uj) = yy{uj) = are zero, then Zy{uj) — deg^(a;), so at both next 
and previous cycle vertex the degree is attained 

deg„„,,,(w) = deg^^^^^Jw) = deg(a;) 

since Zy^^^^ + yy^^^^ ^ Zy + Xy ^ Zy^^^^ + a;„p^^„ (2.38) 

If all degree deficient vertices were of this type, then uj was a multiple of the path consisting of all 
cycle edges. 

Now we divide the set of all deficient vertices (which all have at least one of Xy{uj) or yy{uj) equal to 
zero) into sequences of adjacent ones. Let us fix our attention to such a sequence (we already excluded 
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the case when it has the same end and beginmng). Cah it vi, . . . ,Vr. Then because all w^'s are degree 
deficient: deg(w) — deg^. (w). The last one has yv,X^) > by (|2.38p 

deg(a;) = Zy^_-^{uj) = yvri^) + Zy^i^^) < Zy^i^^) < deg(tj) 

In the same way Xy-^{io) > 0. Finally all the middle ones have Zy.{io) > 1. This implies that we 
can redefine /i on our fixed sequence vi, . . . ,Vr preserving all other properties and fixing the degree 
deficiency: 

2^t)l(A*) = l Xy.{^) = Xy^{fJ,)=0 

Vv, (a*) = yv, in) = yv^ (/i) = 1 

Zyjp)^l ZyX^l)=2 Zy^{/J,)^l 



where i S {2,...,r — 1}. We do this for all such sequences and we have the required /i. 

Now /i is an element of degree 2 and is either a generator or can be decomposed into a sum of two generators 
in degree 1, as described in Lemma [2.371 □ 

Remark 2.39. If we allow more cycles, we can have generators of higher degree. As we can see on Figure [T71 
the graph with two loops and one leaf has a degree three generator: one on the two loops, two on the three 
other edges. 




Figure 17: Example of a degree three generator 



2.5 Embedding. 

The aim of this section is to find a common ambient space for phylogenetic models of all graphs with the 
same topological invariants. The way to construct this ambient space follows easily form the definition of the 
model. 

Theorem 2.40. The phylogenetic model of a trivalent graph Q embeds in a projective toric variety Pg,n; 
which is a good quotient of projective space by an action of a g-dimensional torus. This action as well as the 
variety Pg^n depends only on the first Betti number and number of leaves of Q , up to reordering of coordinates. 

Proof. Models of trivalent trees with N leaves embed naturally in with coordinates where 

K G {0, 1}^ is a sequence of length TV with values in {0, 1} and even number of 1-entries, see Proposition l2.32l 
We cut g cycle edges of the graph Q so that we obtain a trivalent tree T with set S oi M = n + 2g leaves. 
Dividing ~^ by the action of the g-dimcnsional torus T{g, n) that corresponds to gluing these leaves 

back together, yields the required embeddings. 

We will now explain that this action does not depend on the graph Q, up to choice of coordinates on Let 
us label the set of leaves by 1 . . . M and divide it into three disjoint sets 5 = 5*0 U S*! U S"-! as follows. Sq is 
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the set of leaves of the origmal graph Q. The 2g new leaves of T come m pairs {I, l^), where both I and 
used to be the same edge in G- We put I in Si and l~ in S-i. 

The action of the (7-dimensional torus ¥(17, n) on is given by a matrix 

{CiyrelQJ-V ' '^^'^^^ Ck = - Kz- 

Thus this action is independent of the shape of Q and depends only on (5, n) up to choice of order of 
coordinates. □ 

We will now illustrate Theorem 12.401 

Example 2.41. We describe X {^), the model of LittleMan, together with its embedding into Pi, 2- We 
know that the model of a trivalent tree X ^ complete intersection of two quadrics, see |BW07[ Ex. 2.6], 

(xiioQiCooii — a^ooooa^iiii, a^iooia^oiio — a^ioioa^oioi) 

C Proj C[a;oooo5 2;iioo, a^ooiii Xioio, a^iooij 2;oioi: 3:0110] (2.42) 

where each variable corresponds to a degree 1 generator t(X)- By Proposition 12.321 a generator is given 
by its values on the leaves, which form the index k S {0, l}**. 

To obtain X (X), we glue two leaves, say leaf ei is glued with the leaf 62 to get the loop in LittleMan as 
shown on Figure [T51 




Figure 18: Obtaining LittleMan graph from a four-leaf tree 



Now X (X) is a GIT-quotient of X (X) by the action A^^-^^ 

In the same way it embeds in a GIT-quotient of P*". One easily sees that 

PV/C* - V{iji ■y2-zi- Z2) C P(l^ 2^) 

where 

P(l'*, 2"*) = Proj C[a::oooo, a^iioo, a^ooii, aJnn, yi, 2/2, ^1,-^2] 
is a weighted projective space and 

yi — 2:1001 ■ 2^0110, 2/2 — 2'ioio • a::oioi, ^1 = a;ioio ■ a;oiiOj ^2 ~ a;oioi • 2;iooi 

are the C* -invariant variables of degree 2. From Theorem ll.6l it follows that X {^) is given by Equations (|2.42l) 
of X (X) in P^, rewritten in the coordinates of P^/C*. The second one becomes yi — y2 = 0, so 

X(X) = ProjC[T(X)] = ^^(xiiooxooii -2^00003^1111,2/? -^1^2) C F{l\2^) 

II 

n vivi - 2/2) 

n 

P1.2 = PV/C* = Viyiy2~ziZ2) C ¥{l\2^) 
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If we replace each variable in the equations of X(^) by its representation on the graph, we get the picture 
from Figure [TH 




Figure 19: Ideal of the model LittleMan graph 



Example 2.43. As in the previous example, we work out properties of X (o) — the model of Hammock. 
We use the same action A'^~'^, but we change the embedding X Qj^) ^ by relabeling the leaves as on 
Figure HO! 



1\ /3 



Figure 20: Obtaining Hammock graph from a four- leaf tree 

As the labeling of ^ was modified. Equations ()2.42p become 

^(a;ioio2;oioi — 2;oooo2;iiii, a^iooia^oiio ~ ^^iioo^^ooii) C P'^. 
We again rewrite them in the invariant coordinates of P(l^, 2^) to get: 

X{o) = ProjC[T(o)] = V"(a;ooooa;iiiia;iiooa;ooii - ^;i^;2) C F{1^,2'^) 

II 

n ^(yi ^ a^ooooa^iiii, a;iiooa;ooii - 2/2) 

n 

Pi,2 = PV/C* = V{yiy2-ziZ2) C P{l\2^) 

If we replace each variable in the degree 4 equation of X{-o) by its representation on the graph, we get 
the picture shown on Figure 1211 



+ 0+^+^=^+0 

Figure 21: Ideal of the model of the Hammock graph 
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3 Flat families 



Models of trivalent trees that differ by one mutation live in a flat family in a projective space |BW07j . This 
statement almost remains true for trivalent graphs, by the same argument. The only difference is that we 
get a family in the projective toric variety Pg^„ instead of a usual projective space. 

Recall that in Remark 12.161 we have associated to a subset of leaves S C C with k elements an action of 
the /c-dimensional torus T{S) on the model X{0). 



3.1 Key examples. 

We construct two-dimensional flat families containing models of small graphs. They become the building 
blocks for deformations of bigger graphs. 

Example 3.1 (g=0, n=4). Let ^ be a trivalent tree with four leaves. In |BW07| Ex. 2.20] we constructed 
a flat family 

X° ^ Bx P^, 

where 

• S is an open subset of with coordinates 6(i.2)(3.4)j &(i.3)(2.4)7 ^(i.4)(2.3) 

• the torus T(£) acts on B x P"^ via the second coordinate, that is, for a leaf I of ^ and coordinate 
we have A^,(i)(6(.)(.),XK) = t''(')a;K), 

• the equidimentional projection X*^ B contains the three special fibers ^ gj , ^ and A'jq q ^ 
which are models of aforementioned four-leaf trees, 

• X'^ is a T(£)-invariant complete intersection in B x P^ of the two quadrics 

^(i2)(34) • 2;iioo2;ooii + &(i3)(24) ■ 2:^10102:^0101 + ^(i4)(23) ' aJiooi^Jouo 

~ (^(12)(34) + ^(13)(24) + ^(14)(23)) ' a;00002^1111 



(^(13)(24) - &(14)(23)) • a^iiooa^ooii + (&(14)(23) " ^(12)(34)) ' a;ioioa:oioi 

+ (&(12)(34) - ^(13)(24)) • a:;iooia;oiio- 

Example 3.2 = 1, n = 2). We construct a family which contains models of graphs with one cycle and 
two leaves. It arises as a C*-quotient of the family X'^ ^ P'' from Example 13.11 Let us fix a C*-action 
Xh-h ^]-^g ambient P^, thus on X° and tree models as well, by choosing leaves li and I2 labeled by 1 
and 2 respectively. Each of the three trees yields a graph, when two leaves are glued together. Up to graph 
isomorphism, there are two graphs with one cycle and two leaves. As we are mutating along fixed edge, the 
LittleMan appears once and the Hammock twice. The picture of the three possible trees becomes the one on 
Figure!^ 




Figure 22: Mutations of graphs with one cycle and two leaves 
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The new ambient space P''/C* = {yi ■ 1/2 ~ zi ■ Z2) C P(l'', 2"*) was discussed in Example 12.411 

By Theorem ll.6l the new family X'^ //C* is given by the same, T(S')-invariant equations of X^. We rewrite 
them in the (invariant) coordinates of P(l^, 2"^) 

^(12)(34) • 2:11002:0011 + ^(13)(24) • 2/2 + ^(14)(23) ' Hi ~ (^(12)(23) + ^(13)(24) + ^(14)(23)) ' 2;o0002;illl , 
(&(23)(14) - ^(14)(23)) ■ 2;il002;0011 + (^(14)(23) - ^(12)(34)) ' 2/2 + (&(12)(34) " &(23)(14)) ' 2/1, 
2/1 ■ y2 - Zl ■ Z2- 

To understand how this works a little better, let us look at particular coordinate of P(l'*, 2'*), say 1/2 = 
a^iooia^oiiOj and draw on Figure US] its representation for each graph. 




Figure 23: Behavior of coordinates of Pi^2 under mutations 
3.2 General case. 

We construct a flat family containing models of all mutations of a given trivalent graph Q along a fixed edge 
e e £ \ "P. We follow (BW07[ Const. 3.5] replacing polytopes by cones. 

Let Q he a, graph with an inner edge cq which contains two trivalent inner vertices. We can write ^ as a 
sum of a not necessarily connected graph Qi with a set S of k distinguished leaves li,. . . ,lk with fc e {0, . . . , 4} 
and a graph Go having the edge cq as its inner edge and Go is one of the following three small graphs 

(i) a tree ^ with four leaves wi, . . . , W4, where the edge incident to the leaf It is identified with the edge 
incident to the leaf Vi , 

(ii) LittleMan ^ or Hammock -o a graph with four edges and two leaves vi and V2 , 

(iii) Dumbell X a graph with three edges and two loops and no leaves. 

From Proposition l2.25l we can compute the lattice M^^{G) and the cone t{G) from those of the pieces Go and 
Gi 

k 

MS^iG) = (Go) Xg MS^iGi) n fl ker(^, - v,) 

i=l 

k 

t{G) = t{Go) Xg t{Gi) n fl ker(4 - v^). 
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We consider the lattice M^^i^ and the cone Tamb, which are the combinatorial data of the ambient space 
^giGo),n{Qo) described by Theorem l2.401 Example 12.411 and Example 12 .431 The inclusion X{Qo) C ^g{go),n{0o) 
yields maps: M^^f^ M^^{Oo) and Tamb '''(Go)- Forms Vi, i = 1, . . .k pull-back to M^^^^ and we denote 
them by Vi, respectively. Now we define 

fc 

1=1 

k 

Ty = TambiGo) Xg t{Gi) D Q ker(4 - Vi). 



We define a toric variety y = ProjC[rj;]. Since the good quotient is a categorical quotient, by the 
construction we have the embedding 

xiG) y. 

Lemma 3.3. The inclusions 

M§'^M'J^^x Ma-{g,) and Ty ^ Ta^b x riGi) 

induce a rational map 

'^9(oo)MOo) X ^(^i) — y 

which is a good quotient map (of the set over which it is defined) with respect to the action of the k-dimensional 
torus To generated by one-parameter groups Xy.-£., where i = 1, . . .k. The subvariety 

X = X" X X{Gi) ^ BxF'' X XiGi) 

is To-invariant and its quotient X is locally complete intersection in B x y. 

Proof. The map given by inclusions of cones and lattices is a good quotient map by Theorem 11.141 In- 
variance of the resulting subvariety X follows by the invariance of X'~' ^ i3 x discussed in Example 13.11 
Finally, since A* is a complete intersection in S x x X{Gi) its image A" is a locally complete intersection 
in the quotient B x y. This follows from the definition of good quotient, which locally is an affine quotient, 
|BB021 Ch. 5], hence functions defining X locally descend to functions defining X. □ 

Lemma 3.4. Over an open set B' C P^ containing points [1, 0, 0], [0, 1, 0], [0, 0, 1] the projection morphism 
X B' is flat. The fibers over points [1,0,0], [0,1,0], [0,0,1] are reduced and isomorphic to, respectively, 
the geometric model of G and of its elementary mutations along the edge Cq. 

Proof. First we note that the fibers in question, „ oi X B are geometric models as we claim. 
Indeed this follows from the universal properties of good quotients, c.f. |BB02) . as they are quotients of the 
respective products X^^ ^ x X(Gi), which are located, as three invariant subvarieties, m X — X^ x X{Gi). 
This, in particular, implies that the respective fibers of A" — > 6 are of the expected dimension, hence they 
are contained in a set B' C P^ over which the map in question is equidimentional. Since y is toric it is 
Cohen-Macaulay and because A" is a locally complete intersection in 3^, it is Cohen-Macaulay too }Eis95| 
Prop. 18.13]. Finally, the map X ~¥ B' is equidimentional hence it is flat, because B' is smooth, see [Eis95| 
Thm. 18.16] □ 



29 



Theorem 3.5. Geometric models of connected trivalent graphs with the same number of leaves n and cycles 
g are deformation equivalent in the projective toric variety Pg.n, which is a quotient of by a 

g-dimensional torus. 

Proof. This is a combination of Proposition 12.111 and of Lemma 13.41 □ 

Remark 3.6. Since the phylogenetic model of disjoint union of graphs is a product of the models of the 
pieces, see 12.251 we have proved that models of graphs with the same topological invariants are deformation 
equivalent. 

4 Hilbert function of the model. 

As we mentioned in Section|3l given a projective variety embedded in a projective space we have Hilbert func- 
tion coming from the grading of the coordinates ring, or equivalently from the action of the one-dimensional 
torus C*. If our projective variety is equipped with an action of a bigger torus it is natural (see |HS04] ) to 
consider a multigraded Hilbert function, whose domain consists of the characters of the torus. 

We study multigraded Hilbert function of an embedded projective toric variety with the multigrading 
given by a subtorus of the big torus. For the graph model X[Q) the subtorus comes from a subset of leaves. 
We prove in Theorem 14.51 that the Hilbert function depends only on the topological invariants of the graph 
by showing that deformations constructed in the proof of Theorem 13.51 preserve the whole Hilbert series. As 
tools we use Lemma 14.31 to compute the Hilbert series of a torus invariant complete intersection. Lemma 14.41 
gives the formula for the Hilbert series of a quotient of a toric variety by a subtorus of the big torus. We first 
state them in the algebraic setting. 

4.1 Ring with a torus action. 

Let i? be a commutative C-algebra with an action of a torus T. Let Mj = {x : T — > C*} denote the group 
of characters of the torus T. Then we can write 

R = ® Rx 

as a sum of isotypical pieces indexed by the characters of the torus. We assume that each has finite 
dimension over C. Then its Hilbert function iJ/j.T : Mf — > N is 

The Hilbert series is the generating series of Hpi Y 

hRj{t) := ^ dimi?^-t-^. 

XGMt 

Lemma 4.1. Let R = Rx be a ring with a torus action, ///i, . . . , fq are homogeneous with fi G i?^. 

and form a regular sequence in R, then 

hR,Uu...,U)At) = hR{t) .{l-t^-)-...-{l~t^^). 
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Proof. The statement for the single graded Hilbert series is exphcitly given in [Sta78|, Cor. 3.2]. Its 
muhigraded, more general, with minor additional assumption can be found in |MS05[ Claim 13.38]. The 
lemma follows by induction on the length of the regular sequence. For any homogenous / G R^^ which is 
not a zero divisor in R and any x e Af we have the exact sequence of T-modules 

^ i?x-x/ Rx — - Wif))^ — - 0, 

which implies that 

HR/{f)(x) = dim(i?/(/))^ = dimi?x - dimi?^^^/ = Hr{x) - Hr{x - Xf) 

This is equivalent to the required equality for Hilbert series. □ 

Given a subtorus /, : § ^ T we have the corresponding epimorphism of the character groups l* : Ms 
and we can form the S-invariant subring R^ of R equipped with the residual action of the quotient torus T/S 

i?^. R,. 
Then we have the following formula for the Hilbert series. 

Lemma 4.2. Let R be a ring equipped with an action of a torus T and let i : § ^ T be a subtorus. Then 
the Hilbert series of the invariant ring R^ is 

hRSj/s= ^ t^dimR^, 
where i* : Mj Ms is the dual map of the character groups. 

4.2 Toric variety with a distinguished subtorus. 

Let X = Proj i? be a projective toric variety of dimension d with an ample Weil divisor C where R = 
®mGN Rm — ®mGN i "^'^) in Section [Ol Then R has an action of a d + 1 dimensional torus which 

is the product of the d-dimentional torus T of X and the C* from the grading. Any subtorus § ^ T of 
dimension r induces a Z''+^-sub-grading. Then its multigraded Hilbert function Hx,s ■ Mc*xS N 
with respect to S is 

Hx.six) ■■= Hr0{x) = dimi?x- 
The generating series of Hr^s is the multigraded Hilbert series with respect to § 

hx,sit) ■■= hR{t) = dimi?x • t^. 

xeMs 

We have the following corollary of Lemma [4Tl 

Corollary 4.3. Let Y = Proj R be a projective toric variety with an action of an r- dimensional subtorus 
§ C T o/ the big torus. Let us assume that X GY is a ^-invariant complete intersection in Y given by the 
ideal I{X) = (/i, . . . , fq), where deg fi = Xi- Then the S-graded Hilbert series of X is 

hx.s(t)^hY{t)-{l~t^')-...-{l-t^-). 
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The next statement is a corollary of Lemma l4.2l bv using the description of the quotient as the spectrum 
of invariants given in Theorem 11.61 



Lemma 4.4. Let X — Proj R be a projective toric variety with a subtorus t : § ^ T 0/ the big torus as before. 
We assume that both actions are linearized with respect to the ample Weil divisor C Let i* : My — » Afg be 
the correspoding surjection of character lattices. Then 



E 



t^ ■ dim R^ 



Xe-r(X)nkort* 



We combine the above facts to get the equality of the Hilbert series of models of mutation equivalent 
graphs. 

Theorem 4.5. Let Qi and Q2 be mutation- equivalent graphs and Si (respectively S2) be a subset of leaves of 
Qi (respectively S2). Assume that \Si \ ~ \S2\- Then the multigraded series with respect to the tori associated 
to those sets of leaves are equal 



^X(ei),T(5i) 



Jx(e2),T(S2) 



Proof. Since they are mutation-equivalent by Lemma 12.111 we can assume the sequence of mutation takes 
the set Si to the set 32- We can assume Qi and Q2 differ by one mutation. In Section l3^ we have constructed 
a flat family which is a complete intersection having the models X{Qi) and X{Q2) as fibers. Because both 
those models are complete intersections of the same type in the same ambient space by using Lemma 14.31 we 
conclude that the Hilbert series are equal. □ 



We illustrate Theorem 14.51 on examples. 

Example 4.6. We compute the Hilbert series for the models of graphs with two leaves and the first Betti 
munber one X{^) and X(-o-). As we saw in Example 12.341 and in the notation of Example 12.411 the cone 
has generators of the following multidegrees 



coordinate 


2^0000 




2^0011 


3:^1111 


yi 


Zl 


Z2 


degree 


(1,0,0) 


(1,0,0) 


(1,1,1) 


(1,1,1) 


(2,1,1) 


(2,2,0) 


(2,0,2) 



and that X{^) is a complete intersection in P(l^,2'^) of a quadric a;iioo2;ooii ~ 2;oooo2;iiii of multidegree 
(2, 1, 1) and a quartic yf — Z1Z2 of multidegree (4, 2, 2). Thus 

S5) • (1 - t^sfsl) 



^Xa),£a)(^,Sl,S2) - — — .2-^^-, II ,2 2^ n +2 2^ 

(1 - ty ■ (1 - tSiS2Y ■ (x^^-i^iS2) ■ (1 - t^sf) ■ (1 - t^S^) 
Ignoring the multigrading by the two dimensional torus spanned by leaves (by setting si = S2 = 1) we get: 

(1 . (1 _t4) _ l+t2 



hx{X)i't) = 



(1 -t)4 . (1 -t2)3 (l-i)4.(l_t2) 



We have seen in Example 12.431 that t(o) has six generators: those are the elements in the above table for 
^ apart from yi. The model X{-<y-) is a hypersurface of degree (4, 2, 2) in P(l^, 2^) so 



hx{-0),C{-0) {t, Si, S2) 



(1 



(I_t)2.(l_t5^52)2-(l-t25?)-(l-t2s2) 



hxmxm (i, si,S2) 
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Again we can ignore the multigrading and get 



(l-t**) l + t^ 

hxi-<»it) = (I_^)4.(i_i2)2 = (l_i)4.(l_t2)- 

We expand to see the first few terms 

hxa)ii) 1 + 4i + I2t^ + 28t^ + 57t^ + lOAt^ + 176f + 280f + 0(t^). 

Example 4.7. The Hilbert series of models of both graphs and X with no leaves (thus no additional grading) 
and two cycles is 

hxQ)it) = h^^e^it) = _ 4^3 + 6^2 _ 4^ + 1) = 1 + 4t + + 20t' + 35^4 + 56t' 

+ SAt'^ + I20f + 0{t^). 

This is because X{Q) is = (P^ x V^)//{C* x C* x C*). 

4.3 Computing the Hilbert function. 

Given a trivalent tree T with n leaves we computed the Hilbert function Hx(t) of its model in |BW07) as 

where the additional grading corresponds to a distinguished leaf Z, * is an appropriate summing formula and 
Im is the constant function. This inductive formula for Hx{T) uses the decomposition of the tree T as a sum 
of tripods, which leads to the presentation of the polytope A(T) as a fiber product of tetrahedrons A(A,). 

The same method works for any trivalent graph. We proved in Theorem 14.51 that the Hilbert function of 
mutation-equivalent graphs are equal. By Lemma 12.111 we know that any graph is mutation-equivalent to a 
graph of the shape depicted on Figure [24l 




Figure 24: Caterpillar graph 



This means we have reduced the calculation to this case of caterpillar graphs. As we have described in 
Section [2Tl any graph is presented as union of tripods with identifications. More precisely, any trivalent 
graph is built from by the operations of grafting two graphs and gluing two leaves. 

Remark 4.8. To produce a caterpillar graph Q from |'s (leaves) and ?'s (leaves with loop) using ★ and D, 
we need the second operation only once per graph and only in the case when Q has no leaves. 
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Remark 4.9. On the level of graph models we have 
and 

x{g\i-D)^x{g)//c\ 

where the actions of the tori were described in Section 12.41 



We give a formula for /ijt(ei*e2),s(£(Gi*e2)) ^^'^ ^'^^ ^xig'^o) s(£(c;'b)' ^^^^S the above fact about how the 
model of ★ C/2 is constructed from smaller pieces. 



the part containing monomials of the form {tit2t3y{s')'^{s")'~'{s"'y s^oi 

hx(Qt){tl,Si, S„j, — ) ■ h^^J^^{t3,s' ,s" ,s"') ■ hx(g2)it2,Sni + l, ■ ■ .,3,11+712, (4.10) 

where s = (si, . . . , s„2) and / is the exponent vector. 

Let us compute the input functions: apart from the constant one which corresponds to leaves of Q we 
have Hx{i)^s{{i}) the Hilbert function of the model of graph with two edges. Recall that the model X{Y) 
is P^/C* = P(l, 1,2), where the C* action has weights [010 — 1] on P^. Here is the list of generators with 
weights and the resulting graded Hilbert function 







u 


V 


t 


1 


1 


2 


s 








2 



We can expand Formula (|4.10p , setting / to be a Hilbert function of some graph, to get for fc < ^ 



hxi^)*f{k)^{m-k + l) J2 /(a)[2|fc + a](a + l) 



and for fc > ^ 



ao=0 
k 

(fc + 1) f{a)[2\k + a]{m + 2- a) 

aQ—rn — k 
m 

k X! f{a)[2\k + a]{m+l-a) 

aQ—k-\-l 



rn — k—1 



h^* f{k) = (m- k + l) f{a)[2\k + a]{a+l)- 



ao=0 
k 

-fc + 1) Y /(a)[2|fc + a](2m-2fc-a+l) 

aQ—rn — k 
m 

{m~k) Y /(a)[2|fc + a](2TO-2fc-a) 
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where a — (ap, ai, . . . , a„). 

In the same way we can write 



1 



hy,„ii .{t,si,...,Sn) — the part that contains monomials i'(s')*'s'^ of /ix(c') Si, s„, s', —) (4.11) 
where s = (si, . . . , s„) and / is the exponent vector. 
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